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Abstract 

We report a next-to-leading-order (NLO) chiral perturbation theory calculation of 
the neutron-proton scattering cross section in the ^Sq channel using a cut-off regularization. 
The inclusion of two-pion exchanges in the irreducible diagrams - or potential - figuring 
at NLO is found to be important in enlarging the domain of validity of the effective field 
theory. We are able to reproduce the scattering phase shift data up to p = 300 MeV, with 
an agreement which is superior to results of other effective field theory approaches. We also 
discuss the importance of the explicit pion degree of freedom in scattering process. 
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In a series of recent publications, two of the authors (DPM and TSP) together with 
Kubodera and Rho |2|, ^ presented quite successful examples of effective field theory 
in nuclear physics, by showing response functions to electroweak processes at low energy 
involving two-nucleon systems can be described quite accurately with little cutoff depen- 
dence. (For others' works on this area, see Refs. [@, ||, |6|, |^, ^.) The approach used there was 
the "A counting scheme" proposed by Weinberg The focus on the response functions 
instead of on scattering observables is in the spirit of the long-standing tradition in nu- 
clear physics where the wealth of nuclear dynamics has been accumulated more powerfully 
through response functions than through scattering processes per se. 



The recent surge of activity in effective field theories in nuclear physics |1C] was 
however triggered by the effort to give a first-principle description of the large scattering 
lengths in nucleon-nucleon scattering, so a large portion of the recent publications has 
been devoted to scattering amplitudes at low energy. "Q counting" scheme introduced by 



Kaplan, Savage and Wise |11| using power divergence subtraction (PDS) exemplifies the 
present preoccupation of the workers in this field. ^ While there appeared a few papers 



on the power of the A scheme in scattering p3| , 14, 15], an exhaustive confrontation of the 
A scheme developed for the response functions Q with scattering amplitudes has not been 
performed. It is the purpose of this paper to provide the missing information. We shall 
show that the scheme is as successful for scattering as for response functions and in fact is 
even more accurate than the Q counting scheme. 

The question we address is this: How does the A scheme of Q which is remarkably 
successful in postdicting and predicting electroweak processes of two-nucleon systems p+p —>■ 
d + e'^ + V [Q] and n + p— >(i + 7 [^, fare with the two-nucleon scattering and to what 
nucleon momentum and with what accuracy can one "push" the scheme? 

For np scattering to the next-to-leading order (NLO) in the A scheme with pions and 
nucleons, the two-particle irreducible graphs that generate the potential for the system with 
which the Schrodinger(or Lippman-Schwinger) equation is to be solved comprise of nucleonic 
contact interactions and pion-exchange potentials. The contact terms, the coefficients of 
which will be fit to empirical data at the zero momentum limit, represent the short-range 
part of the nuclear interaction. On the other hand, the pion exchange potentials control 
the long range part of the nuclear interaction, so unsurprisingly the explicit inclusion of 
the pion-exchange potential generally increases the range over which the wave functions are 
accurate. It was shown in that the absence of the OPEP could barely be compensated 
by the leading-order contact term and that the contact interaction at the next-to-leading 
order (NLO) - when added to the LO interaction - is seen to improve considerably the 
stability of the numerical results, thereby increasing the domain of validity of the effective 
field theory. In this paper, we will take explicitly into consideration the two-pion-exchange 
potential(TPEP) at the NLO and show how accurately one can calculate the two-nucleon 

*^ There is also a slightly different scheme which emphasizes the unitarity and the relativistic formalism 
suggested by Lutz |l2[. 
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Figure 1: Two-body irreducible diagrams that contribute to the TPE potential. 



process. In addition, we will also report on a higher order calculation. 

We have described in great detail the main strategy of the cut-off EFT in our previous 
works, so we will not enter into it here. It will suffice to briefly define the convention and 
go directly into the results. We adopt the power counting rule given in ||9|: an irreducible 
diagram is of order of 0(XQ I ^xt") where Q is the pion mass or the typical momentum scale 
of the process and ~ 1 GeV. And we keep only pions and nucleons as pertinent degrees 
of freedom, all other massive degrees of freedom are integrated out. In nucleon-nucleon 
potential, the LO has z/ = and consists of OPEP plus non-derivative contact interactions: 

VLo{q) = TjCo + ViAo), Vi.(g-) = 1^^^, (1) 



where M is the nucleon mass, — 1-25 the axial-vector coupling, ~ 93 MeV the pion 
decay constant and q the momentum transfer. There is no = 1 contribution, the NLO 
(v = 2) comes from TPEP and contact interactions with two derivatives. TPEP contains 
also v = 3 (and higher order) contributions, which correspond to the subleading order 



ttttNN vertices. There have been some elaborate works on the TPEP [17, 18, 1£] up to 
v = 3 order and applications of the potential in other channels as well as the S-wave np 
scattering. In this work, we will not consider the = 3 order contributions, though it may 
results in substantial improvement of the theoretical predictions. 

The two-nucleon irreducible diagrams that contribute to the TPEP are shown in 
Fig. In addition, we also include the contributions from OPEP-subtracted two-pion-box 
diagrams, whose effects are negligible in the ^5*0 channel. Remembering that the contact 
interaction with the coefficient C2 enters at the same order, we write the full NLO potential 
as 



with 



647r2/4 



dx 



4m| 



+ 



f2n{x) 



4x^ + 8ff^x^(3 + 2x^) - 4ff^(3 + 48x^ + 8x'^) 
24(1 - x2) ' 



(2) 



(3) 
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To the same order in the chiral counting, there are loop corrections in the propagators 
and the ttNN vertex. However in the kinematics involved in the elastic scattering, they 
do not bring in any momentum-dependent corrections: They modify only the constants 
(masses and coupling constants) which are to be absorbed into the parameters extracted 
from experiments. |^ 

So far, we have assumed that \q\ is of the same size of the pion mass m.,^. When 
\q\ is much smaller than myr, we can integrate out even the pions. Theory then contains 
only nucleons explicitly, and the expansion parameter becomes \q\/m.,^, rather than Q/A^.. 
As we see, OPEP and TPEP become of order of |g|^/m^ and |g|^/(2m^)^, so they can be 
replaced by suitable series of contact interactions starting from two-derivatives and four- 
derivatives. By the same token, multi-pion-exchange potentials can also be replaced by 
contact interactions. Recalling that multi-pion-exchange potentials are hard to calculate, 
this procedure provides us an efficient way of going to higher order. The price of this 
simplification is the shrinking-down of the radius of convergence. Following to the general 
argument of the effective-range expansion, the radius of convergence is about a half of the 
lightest mass of the degree of freedom that is not explicitly taken. That is, without the pion 
degrees of freedom, we can describe only up to p < where p is the CM momentum. 

Taking OPEP explicitly with higher order contact interactions, we then expect that we 
can achieve accurate description of the scattering data up to p < m.,^, which was indeed 
proved in |||]. In this paper we go to one further step, by exploring the consequences of the 
inclusion of the TPEP. With the TPEP, we expect we can go up to p < ^m-n^. Basically, in 
this work, we will focus on the comparison between the NLO potentials with and without 
the TPEP. For clarity, we shall refer to the latter as the 1ttC2, meaning that it contains 
OPEP and contact interactions up to C2 term. Thus in this notation, LO is identical to 
IttCq and NLO is identical to 27rC2. One of our principal findings from this comparison 
between the NLO and l7rC2 potentials is that the explicit presence of the pion within the 
A scheme makes the EFT a lot more versatile and accurate. In fact, this work provides 
a confirmation of the partial finding in Q of the effect of the explicit OPEP in NLO and 
1ttC2 in scattering consideration. We will also consider 27rC4 potential, which consists of 
the NLO potential and four-derivative contact interaction of the form ^C4|g|^. This 27rC4 
potential as a possible next subleading order may only find its justification for \q\ smaller 
than m,r) whereas for other kinematical zone, the inclusion of all subleading NNLO TPEP 



should be made |21|. Thus compared to the NLO, we expect that the theory can improve 
(both in accuracy and in cutoff-independence) results at that low momentum region, but 
not the radius of convergence in general. These hypotheses will be proved in our work. For 
the low momentum region with p <^ 1^,^, we expect this potential provides us a simple 
but sufficiently adequate way to do an 0(|q|^)-order calculation without including pion 
two-loops. It is worth noting that the integral in Eq.(^) can be rewritten in a spectral 

*^ There are also certain contributions that come from relativistic corrections. A detailed analysis of these 
contributions is recently given by Friar |2d]. Here we are not concerned with them. 
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representation by changing the integration variable j; to a "spectral mass" m ~ ^™ 



which runs over 2m.j^ < m < oo. One may then introduce an upper limit to the spectral 
mass m < uiq, by imposing the range of the integration variable to be < x < wl — 
The mo should not be confused with the genuine cutoff of the theory A that enters through 
the Fourier transformation as defined below. Note that the TPEP is negative while OPEP is 
positive for the whole range of \q\ and that it is quite small compared to the OPEP, less than 
3 % for the momentum range \q\ < m-n. This implies that the results obtained in will not 
be changed appreciably by the explicit inclusion of the TPEP. As the momentum increases, 
however, the genuine TPEP will play an important role for the radius of the convergence, 
the cutoff-dependence and the allowed range of the cutoff of the theory. Throughout this 
work, we put the value of mg large enough, mg = 2 GeV - instead of infinity - to economize 
the calculation. 

These potentials of various chiral orders, transformed into coordinate space, will be 
put into the Schrodinger equation. In doing the Fourier transform, a cutoff in momentum 
space will be introduced. The potential with the cutoff A takes the form 

V{f^ ^ / (0 ^^'"^ ^A(g') V{q) (5) 

with 

S^{<f) = e-^. (6) 

Due to the cutoff - which is kept finite, the integrals are finite so that we don't need counter 
terms in a proper sense in addition to the contact terms that figure to the order considered. 
However, the concept of renormalization is still valid and figures in the requirement that 
the results be stable against the variation of the cutoff. The Fourier transformation of the 
OPEP and the contact terms can be done straightforwardly (see, for example, jl], ^). For 
the TPEP, a similar technique can be used by performing the Fourier transformation before 
the parametric space integration. The wave function ^o{r) = is then obtained through 
the Schrodinger equation 



^^+M{E-V{r)) 



uoir) = (7) 



where E = /M is the total energy of the system, and p = \p\ is the center-of-mass 
momentum. The solutions are to fit the empirical effective range expansion, 

1 1 

pcoi5 = V -reP^ + V2P^ + 0{p^). (8) 

«o 2 

For this we shall take the Nijmegen partial wave analysis (NPWA) [^] as the empirical 
data and extract from it the values of the low-energy constants, = —27.73 fm and r^. = 
2.677 fm. Note that the NPWA results are in very good agreement with the values obtained 
by the Argonne fis potential [p^j , i.e., ao(wi8) = —27.732 fm and re(fi8) = 2.697 fm. The 
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0.14 


0.38 


1.56 







Table 1: Coefficients of the contact interaction terms in LO, lvrC2, NLO and 27rC4 poten- 
tials, for different values of the cutoff A. We have multiplied A" to the Cn {n = 0, 2, 4) so 
that the unit of the Cn-^"" is fm for any n. The blanks for the 27rC4 are explained in the 
text. 



constant C4 will be determined to reproduce the effective volume V2 = —0.48 fm^ |2^, which 
was extracted from the NPWA data [23|.^ We are now in position to discuss our results. 

In confronting the experiment with our theory, we shall particularly be interested in 
assessing how well the basic tenets of an effective field theory of the type we are considering 
are satisfied |22|. The first is that the effective field theory in question is a nonrenormalizable 
theory and a truncated version with a given cutoff scale must break down at some point 
signaling the emergence of "new physics." The second point is that the effective theory with 
a cutoff regularization is valid only if the physical observables calculated within the scheme 
become stable against the variation of the cutoff when higher order terms are included, which 
is nothing but an approximate manifestation of the renormalization group (RG) invariance 
that any valid effective theory must satisfy. 

Let us see how these points emerge in the calculation. In Table H are displayed the 
values of the coefficients for various potentials. We can observe in the Table that, except 
for the LO case, the values of C's increase rapidly at a certain value of the cutoff. In fact 
there is a critical value of the cutoff beyond which we cannot find the Cs that reproduces 
the low-energy constants |26|. The 27rC4 for A > 500 MeV is found to be such a case.^ We 



show in Fig. ^ the ^Sq phase shift as a function of the cutoff A for the initial CM momenta 

of 70, 140, 210 and 280 MeV (from top to bottom) for different potentials. In the figure, 

the NPWA data are represented by thin horizontal solid lines; the LO, 1itC2, NLO and 

*^In jl^, a slightly different renormalization procedure has been introduced. They are using the modified 
effective range expansion, which deals with the OPEP-subtracted A'^A'' potential. Their method is not 
adopted in our work. 

If we allow a slight deviation by about 0.03 fm"* in the effective range, we can go at least up to A = 600 
MeV. 
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27rC4 results by dotted, dashed, solid and dotted-dashed curves, respectively. The extensive 
comparison between LO (dotted curves) and 1ttC2 (dashed curves) in this framework has 
been given in [|^: By adding C2 term into the LO potential, the l7rC2 theory becomes 
much more accurate with much less cutoff dependence compared to the LO theory. Indeed, 
even at p = 140 MeV where the theory is expected to break down, the theory differs from 
experiment only by less than 3 % for all the region with A > 300 MeV. We also observe 
that the theory becomes closest to the experiment at A ~ 300 MeV, while the LO favors 
A ~ 200 MeV. Now by adding the TPEP, we see that the NLO (solid curves) becomes more 
accurate than the 1ttC2 at higher cutoff. For all the region considered with A > 300 MeV, 
the NLO differs from experiment less than 13 % at p = 210 MeV. Recalling that the 
theory is expected to break down for p > Im^r, the above result is encouraging. It is also 
worthwhile observing that the NLO with A = 400 MeV nearly coincides with the NPWA for 
all the momentum considered, p < 280 MeV. On the other hand the improvement between 
the l7rC2 and the NLO is not spectacular. The fact that the 27rC4 cannot come closer to 
the NPWA shows the importance of the omitted subleading TPEP and higher order terms. 
Thus mere inclusion of higher-order contact interactions cannot possibly ameliorate the 
convergence. This indicates that any attempt to improve the precision by merely including 
higher order contact terms without the explicit account of the relevant degrees of freedom 
that are nearby is doomed to fail. 

In the figure, it is clearly seen that the theories with higher terms become more 
insensitive to the cutoff, which shows the realization of the above mentioned RG invariance. 
At p = 280 MeV, while the NLO shows reasonable agreement with experiment, results show 
large fluctuations with respect to cutoff which implicates the breakdown of the theory. From 
Fig. 3, one can confirm that the theory breaks down at around this momentum. 



For completeness, we shall now compare our cutoff scheme with the PDS scheme[27, 
28 1 . For ease of comparison, we estimate the deviation from the NPWA result for the 
effective range function cot 6 of the cutoff EFT and PDS schemes as a function of the 
momentum. The results are shown in Fig. ^ Our NLO results are drawn in solid curve 
and the l7rC2 in dashed curve. The PDS results obtained up to NLO are drawn both for 
the coefficients fit to the low-energy constants, oq and r^. (dotted line) |l^ and for the 
results obtained by a global fit for the range p < 200 MeV (dot-dashed line). It might be 
worthy to note that the NLO in PDS and in our framework have different meanings, as the 
counting rules adopted are different. Nevertheless both theories have the same number of 
free parameters, so the comparison can be meaningful. From the figure, one can confirm that 
the cutoff theory is doing better than the PDS results in accuracy. A similar observation 
was made in |14]. Recently, the NNLO calculation in the PDS method (which is order 
compared to the LO) was reported, both in a toy model |27] and in a real situation |29|. 
In [^], they have a good agreement up to p = 400 MeV by fitting the parameters over 
p = 7 — 200 MeV. Due to their fitting procedure, it is however meaningless to directly 
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Figure 2: Phase shift with respect to A at p = 70, 140, 210 and 280 MeV. For each 
momentum, the LO, l7rC2, NLO and 27rC4 results are represented by dotted, dashed, sohd 
and dotted-dashed curves. 
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Figure 3: | A cot 5\ (= | cot (5epx ~ '^NPWA I) with respect to the momentum p. SoHd line 
is for NLO(= 27rC2) and dashed hne for l7rC2. Dotted line is the results of PDS fitted to 
oq and Tg and dot-dashed line is the results for fit over p < 200 MeV. 



compare their results to ours.^ In applying their framework to the spin-one channels, they 
conclude that the nonperturbative treatment of potential is needed. 

To conclude, we have demonstrated in this letter that the A scheme can be made 
markedly more successful for the scattering problem by going to the NLO order a la Wein- 
berg. We have also shown the importance of the explicit presence of the pion exchange 
potential in NLO, extending the domain of validity as well as increasing the precision. The 
method is found to be as successful in scattering as in electroweak responses. It has also 
been discussed that there are missing parts in the TPEP, the subleading order contributions, 
which are expected to play a significant role in increasing the accuracy at high momentum 
region. While we are revising this paper, quite a similar work has been reported by Epel- 
baum, Glockle and Meifiner |21]. They have constructed the potential including up to the 
subleading TPEP (and the leading TPEP but with also Delta isobar). They have performed 
the calculation in momentum space, and applied to various channels as well as the l^o np 
scattering. And they have got similar conclusions to ours, including the importance of the 
subleading TPEP. 

We are particularly grateful to Mannque Rho for his continuous supports and invalu- 
able discussions during this work. Work of DPM and CHH is partially supported by the 
KOSEF through OTP of SNU and by the Korea Ministry of Education under contract No. 



*^ As a consequences of the fitting, their low-energy constants are different from experiment. For example, 
the effective volume in their approach is V2 = —1.2 fm^, which should be compared to the data —0.48 fm^. 
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